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Abstract 

We study transport equations for quantum many-particle systems in terms 
of correlations by applying the general formalism developed in an earlier pa- 
per to exactly soluble electron-phonon models. The one-dimensional models 
considered are the polaron model with a linear energy dispersion for the elec- 
trons and a finite number of electrons and the same model including a Fermi 
sea (Tomonaga-Luttinger model). The inclusion of two-particle correlations 
shows a significant and systematic improvement in comparison with the usual 
non-Markovian equations in Born approximation. For example, the improved 
equations take into account the renormalization of the propagation by the 
self-energies to second order in the coupling. 

I. INTRODUCTION 

In a recent work |]J one of the authors proposed a new method for deriving transport 
equations in terms of correlations. In this paper we test these transport equations and 
gain more insight into their abilities to describe physical effects for some exactly soluble 
electron-phonon models. 

In their study of the Jaynes-Cummings model Zimmermann and Wauer used the 
same decoupling procedure as ours one order higher than the Born approximation in order 
to reduce unphysical results such as fermionic occupation numbers outside the range [0,1]. 
We extended this procedure systematically to all orders in the correlations, thus obtaining 
an exact, infinite set of equations of motion in terms of correlations. For this simple model 
it is easy to solve these equations numerically to very high order in order to test our method 

Here, we focus on one-dimensional versions of more realistic models, i.e. the polaron 
model and the same model with a Fermi sea of electrons (Tomonaga-Luttinger model). 
For the polaron model with a linearized electronic energy dispersion exact results for non- 
equilibrium expectation values were presented in Ref. H, as well as for the Tomonaga- 
Luttinger model with coupling to phonons in Ref. ||. Similar models have recently been 
discussed to study the relaxation of highly excited semiconductors and insulators 
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For these models non-Markovian quantum-kinetic equations in Born approximation were 
numerically solved assuming isotropic distributions in fc-space. It became clear that damping 
factors for the electron propagation had to be introduced in order to reduce the occurrence 
of unphysical results, e.g. negative occupation numbers. Renormalized free propagators 
including the damping were obtained by the generalized Kadanoff-Baym ansatz fllUyTyTTf , 
which in turn was justified by Schoeller |12| as a partial resummation in a diagrammatic 
perturbation theory. For the Tomonaga-Luttinger model including phonons it turned out 
that the renormalization of the phonon frequencies plays a major role for the relaxation 
process ||. In this work we study in which way these important renormalizations of the 
free propagators appear in our method, which contains only single-time quantities without 
using advanced and retarded Green's functions. This is quite interesting as the Hamiltonian 
dynamics is local in time and memory effects only appear by truncating the infinite hierarchy 
of quantities. 



II. THE MODELS 

As in the first paper, we consider systems described by a general statistical operator p 
at the initial time t^. The dynamics is determined by a Hamiltonian H, and the expectation 
value of an operator A at time t is given by {A) t = Tr[p A(t)], where A(t) is the operator 
in the Heisenberg representation with A(to) = A. 

The one-dimensional polaron model is described by the following Hamiltonian 

h = E £k4i>k + E ^ B l B i + £ & (4 +q ipkB q + B\4^ k+q ) . (i) 



k q q,k 

The operator ip k (B q ) annihilates an electron (phonon) in the state with momentum k (q) 
and energy e k {yj q ). The electron-phonon coupling strength is denoted by g q . In the case of 
a linear energy dispersion e k = vpk for the electrons this model is exactly soluble for the 
initial one-electron state po = l^oX^ol or more generally for p = • • • , fcjv)(&i, • • • , &at| 
where \ki, . . . , k^) is an anti-symmetrized iV-electron state QT3 1 . The time evolution of the 



electron distribution function n k (t) = (V'lV'fc)* for the initial state p = \k )(k \ can, e.g., be 
obtained by numerically solving the exact differential equation derived in Ref. ||: 

d , s ^ o sin \{uj a — VFQ){t — t )l , . 

-n k (t) = 2j29 2 q [V \ ™ °-^[n k+q (t)-n k (t)] . (2) 

at q>0 u q — VFq 

As in Ref. Q we assume that only phonons with q > exist in accordance with the fact 
that each electron is moving "right" with the velocity Vp (chiral model). Even for the initial 
state po = \ko] FS)(k ; FS] (k > k F ) with a Fermi sea |FS) = J2k<k F ^IlVac) of electrons it 
turns out that the exact time evolution of the electron and phonon distribution functions 
can be numerically determined |J. Here we use the exact solutions for a comparison with 
the results of our approximate transport equations. 

The exact solutions are obtained for the case of an infinite range of momenta k = 2rm / L 
with n = . . . , —1, 0, 1, . . . whereas the approximate numerical calculations are for a system 
with a lower and upper momentum cut-off (k < < k < fc>). These differences are of little 
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importance if the cut-offs are chosen to be far away from the Fermi momentum resp. the 
initial state at kg. In all numerical calculations we use a constant electron-phonon coupling 
with a momentum cut-off q c , i.e. g q = gQ(q) Q(q c — q), and a constant phonon frequency 
ujq = ujg as a model for longitudinal optical phonons. We define the "resonant" phonon 
momentum q B by u = Vpq B and choose q c ^> q&. In the plots of the numerical results 
all momenta k will be given dimensionless as k/q B and all times t as uiot- The system 
is described by the following dimensionless parameters: the ratio of the one-electron level 
spacing and the phonon energy v = (v f^tt / L) / u = (2ii/ L)/q B , the coupling constant 
r = g 2 /(u ■ v f 2tc/L), and q c /q B . 

III. QUANTUM KINETIC EQUATIONS 

We now present the equations of motion (EOM) for the correlations derived using the 
general formalism described in Ref. |TJ. First, the transport equations one order beyond the 
Born approximation are studied. The form of the equations is independent of po, which only 



determines the initial values. The relevant correlations are (q,q' > 0): 

n k {t) = (4i> k ) t (3) 

N q {t) = (B\B q ) t , (4) 

h k , q (t) = i ([B\4_^ k ) c t - (B q 4^ q )fj (5) 

H k , q (t) = (Bl4_ q ip k ) c t + (B q 4^ q ) c t , (6) 

n k>q (t) = 2(4 k _ q 4x!> k _M c t (7) 

n k , q>q >(t) = (tpl-q-q^Uk-q^k-q'Tt + {ii-q'^k-q^k-q-q'Tt (8) 

N k , q , q >(t) = i ((^i-q-qi^k-q^k-q'Yt ~ {^k-q'^k-q^k-q-q'Yt) , (9) 

*W(*) = (B\B q ,4_ q ^ q ,) c t + (Bl,B q 4_ q/ ^ q ) c t (10) 

K k>q>q ,(t) = i ([B\B q ,4- q ^k-q>) c t - {Bl,B q 4_ ql ^ q ){) , (11) 

l k , q , g/ (t) = [B\B\4\^ k ) c t + {B q ,B q 4^ q _ q ,) c t (12) 

L w (f) = i {{B\Bl,4_ q _ q 4 k ) c t - (B q ,B q 4^ k _ q _ q ,){) . (13) 

The indices are chosen such that the correlations with three indices are either symmetric or 
anti-symmetric in q and q'. The transport equations read 

4: n k(t) = J2dq (h k>q (t) - h k+q , q (t)) , (14) 

at q>0 
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l Nq (t) = -g q J2h k , q (t) , (15) 

k 



j£ h k,q(t) = (Cfc - e k _ q - UJ q )H k;q (t) (16) 

+2^ [jV,(t) (n*_ 9 (f) - n k {t)) - n k (t) (1 - n k - q (t))] 
-9q[n k , q (t) + + n k+q ,^ q ,(t))\ 

q'>0 

+ H fiV (h,q,q'(t) - h+q',q,q'{t) ~ k,q,q'(t) + k+q',q,q> (*)) , 
g'>0 



^H k>q (t) = -{e k - e k - q - u q )h k , q (t) (17) 

~09 S ( N k,qAt) - N k+q',qA t )) 
q'>0 

- 9q' {K k , q , q '{t) - K k+q ,^ q ,(t) - L k ^ q ,{t) + L k+q ,^ q ,(t)) , 

q'>0 



^ n *.«(*) = 2 #7 ( n k- q (t) - n k {t)) h Kq {t) , (18) 

-j: n k,q,q'(t) = {e k - q - q ' + Cfc — — € k - q >) N k q q '(t) (19) 

at 

+ {#9 [(n fc _ g _,#(t) - n k-q'{t)) h k>q (t) + (n k -q(t) - n k (t)) hk-j,q(t)\ 

-(?W)} , 

~^^k,q,q'{t) = —{^k-q-q' + ^k ~ e k-q ~ £k~q') n k,q,q' (t) (20) 

+ [(^fe- 9 -?'(i) - n k - q/ (t)) H k>q (t) - {n k _ q (t) - n k (t)) H k _ q/>q (t)] 
-(?W)} , 

j^k k ^ q/ (t) = (e k _ q + u q - e k _ q i - uj q ,)K k ^ q i{t) (21) 
+ [(JV,(t) + n w (t)) /i fc _^(t) - (1 + W,(t) - ra fc _,(f)) V'Wl 

+(?<-> 5')} , 

jj.Kk,q,q'(t) = —(Ck-q +^q~ ^k-q' ~ Vq')kk,q,q'(t) (22) 

+ [(N q (t) + ra w (i)) - (1 + N q (t) - n k _ q {t)) H Kq ,{t)\ 
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-^h,q,q'(t) — (&q + <*ty + ^k-q-q' ~ e k) Lk,q,q' (t) (23) 

+ {g q [- {N q {t) + n k (t)) h k ^ q> (t) + (1 + N q {t) - n fc _,_g/(t)) h k , q ,(t)} 

~Jj.Lk,q iq '(t) = —(uJq + UJ q i + tk-q-q' — £k)h,q,q'(t) (24) 

+ {g q \(N q (t) + n k (t)) H k „ m ,{t) - (1 + N q (t) - n k - q -g>(t)) 
+{q <-> g')} • 

In the graphical representation of the EOM Jl| we often use fermionic (solid) and bosonic 
(curly) lines without arrows in order to reduce the number of diagrams. Arrows and mo- 
mentum labels are supposed to be added according to the conservation of the number of 
electrons and of the momentum. The interaction vertices are given in Fig. [I]. The graphs for 
the distribution functions n k and N q are those shown in Fig. |2|, the diagrams for the "phonon 
assisted polarizations" h ksq and H k q are shown in Fig. |3], the graph for k k>qiq i, K kjq>q i, lk, q , q ' 
and Lk t q <q ' is given in Fig. f| and the graph for n k>q , n k ^ q t and N k , q , q i is the one shown in 
Fig- B 

By restriction to the correlations n k , N q , h k ^ q and H kq the transport equations in Born 
approximation are obtained. For the case of vanishing initial correlations they read after 
integrating the differential equations for h k q and H k q : 



dt 



n k {t) = 2Y.9q f dt' {cos ((e* - e k - q - u q ){t - t')) (25) 
x [N q (n fc _, - n k ) - n k (1 - n k _ g )} t , 

- cos ((e k+q - e k - LUq)(t - t')) [Nq (n k - n k+q ) - n k+q (1 - n fc )] t ,j , 

^-Nqit) = -2g 2 q Yl f dt' cos {{e k - e k ^ q - u q ){t -t')) (26) 
at k Jto 

x [Nq {n k _ q - n k ) -n k (l- n k _ q )] t , . 



A. The polaron model 

We first discuss the transport equations for the one-electron polaron model with the 
initial state po = |&o)(&o|- A renormalization of the free propagators, i.e. the cosine factors 
in Eqs. (|25|) and (p6|) , is of importance for the Born approximation to yield acceptable results 
for the first phonon replica at larger times [||,[7| . In the one-electron case there is no significant 
renormalization of the phonon propagators and the most important renormalization is the 
inclusion of a damping factor exp[— 2^{t — t')] for the electronic propagations. For the 
determination of 7 from the second order electron self-energy in equilibrium we refer to 
Ref. [§. 



5 



We solved numerically the equations including the two-particle correlations as described 
in Eqs. (|T4l)-(|2~4]). The initial states discussed so far do not have initial correlations. Fig. [7] 
shows the time evolution of n k for the first phonon replica at k = k — q B with q B = u> /vp. 
The importance of the damping of the electronic propagation in the Born approximation can 
be seen clearly. The approximation including the two-particle correlations seems to include 
the damping in the dynamics as it follows the exact evolution for a long time. 

We now explain how the dynamical equations incorporate the damping of the electron 
propagation. In the Keldysh formalism |14| the renormalization of the free propagators is 



usually justified by the Kadanoff-Baym ansatz |TU|-[T^], which substitutes the free propa- 
gators in Eqs. (^) and ([26|) by the exact retarded/advanced Green's functions. As these 
Green's functions are not known they are, in a first approximation, taken in lowest order per- 
turbation theory for the equilibrium. In our case, the electron propagators are renormalized 
by the self-energy in second order Born approximation shown in Fig. || Diagrammatically 
these self-energy insertions arise from the third diagram in Fig. Therefore, in order to 
explain the renormalization we can neglect for a moment the two-electron correlations n k S1 
nk^q^q' and Nk } q, q >. Then the differential equations for the phonon assisted density matrices 
read 

= 2g q [N q (t) (n fc _,(t) - n k {t)) - n k (t) (1 - n fc _,(t))] 

+ ^2 9q' ((kk,q,q' ~ ^k,q,q' ) (t) ~ (kk+q>,q,q> ~ iKk+q' ,q,q'){t) 
q'>0 

— (lk,q,q' — ^k,q,q'){t) + Q-k+q' ,q,q' ~ ^fc+«' )«>«') (*)) 

The terms in the time derivative of (k ^ — iK ) and (7... —iL ) that yield the renormalization 
of the electron propagation for (h k}q +iH k}q ) are the "coherent terms", i.e. those proportional 
to h k , q and H k>q : 

~~ ^ Cfc -<? + UJ 1~ €k -l' ~ ^V)^ ( k k,q,q> - iKk,q,q'){t) (28) 

= -g q , (1 + N q ,{t) - n k - q ,(t)) {h hq + iH ktq )(t) +... 



-Jl - i( e k+q'-q + - £k - ^g')^ (h+q' ,q,q> ~ ^fc+g'.g,?') (*) ( 2 9) 



At 



= g q < (N q ,(t) + n k+ g'- q (t)) (h Kq + iH k>q )(t) + 

and the analogous expressions for (I —iL ). If there are no initial correlations we obtain the 
following integro-differential equation for (h k>q + iH k ^ q ) neglecting the "incoherent terms": 

2g q [N q (t) (n k - q (t) - n k (t)) - n k (t) (1 - n k - q (t))] 



f ' dt'U k , q (t,t')(h k , q + iH ktq )(t') 

Jto 
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where 

U kiq (t, = Q(t - 2 <7 ? 2 , [(1 + N q , - n w ) {t'y^-^-^'-^-^ (31) 

q'>0 

+ (N q/ +n k+q/ ) (t') e <^+^+^- g -e k+q r)(t-t') 

+ (1 + N q , - n fe _ g _ 9 (t') e i K+^'+^- 5 - 9 '-^)(*-*') 

+ + n fe+g ,_ 9 ) (^(w-^-^-vX*-*')" . 

In this context, the time- dependence of the distribution functions n k and N q is of little 
importance, as can be checked numerically. If we therefore replace them by their initial 
values rife(to) an d N q (t ), the "self-energy" H k)q (t,t') depends only on the relative time t — t'. 
Then the differential equation (^) is easily solved: 

(h kjq + iH k , q )(t) = 2g q f dt' D k<q (t - [N q (t') (n k - q (t f ) - n k (t')) - n k (t') (1 - n k ^ q (t'))} 

(32) 

where D k , q is the retarded solution to the Dyson equation 

D Kq (t - t') = Dfl(t - t') - £ dt x £ dt 2 D^(t - ti)n M (t! - t 2 )D k , q {t 2 - f) (33) 

with the "free propagation" D^(t - t') = 0(t - t')e-^ e * - e *-«" w «) (*-*'). Except for the fact 
that D k q is the "propagator" for an electron-hole pair plus a phonon, this is an ordinary 
Dyson equation, which can be easily solved by Fourier transformation 

- iD kA (u) = (u - (ejfe - e fc _ 9 - u q ) + iU k , q { u ) + ( 34 ) 

and leads to a damping for both electronic propagations. 

As can be seen from Fig. [8] the time evolution of the occupation at the second phonon 
replica at k = k — 2qs is much better described by the approximate equations including 
correlations (long dashed line). This is easily explained by the fact that contrary to the Born 
approximation these equations contain coherent two-phonon processes by the inclusion of 
expectation values [B^'B^'^ip)^. By an analogous argument to the one used to explain a 
damping factor in the Born approximation it can be justified to damp the time evolution 
of n k , q , n k}m >, N k>qjq >, kk, q , q >, Kk,q,q>, h,q,q', and L k , q , q >, i.e. replacing j- t by (jg + 27) in the 
corresponding differential equations. Why the correlations with four electron operators 
(ijj^ijj^ijjijj)^ are damped with 27 instead of 47 as one might expect will be discussed later 
when we extend the equations by "half an order". The inclusion of this damping for the 
two-particle correlations yields an improvement of the solution for longer times in the same 
way as the damping of the Born approximation does one order higher (dot-dashed line in 
Fig-i- 

Figs. |9| and |10| show the electron distribution function at two different times. In the 
range of the first phonon replica, the distribution function is well described by the Born 
approximation with a damping factor for the free electron propagation (dashed lines) and 
by the kinetic equations including two-particle correlations (long dashed lines), whereas the 
Born approximation without damping (dotted lines) yields unphysical results like negative 
occupation numbers. For the second phonon replica only the kinetic equations including 
two-particle correlations are in excellent agreement with the exact solution (solid lines). 
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B. The Tomonaga-Luttinger model 



We now turn to the model with a Fermi sea (Tomonaga-Luttinger model). The initial 
state p — \ko] FS)(k ; FS| again contains no correlations. The most important effect of the 
Fermi sea (apart form the Pauli blocking) is the renormalization of the phonon propagator. 
Due to the linear energy dispersion of the electrons, which also enables the exact solution 
by bosonization || , the phonon mode splits into two modes in the equilibrium case and the 
second order self-energy (electron-hole bubble) already gives the exact result for the equi- 
librium phonon Green's function. In Ref. || the importance of the phonon renormalization 
for the kinetic equations in Born approximation was shown. In our dynamical equations 
this renormalization is incorporated by the second diagram in Fig. |3|. By extending the dis- 
cussion following Eq. ^ to the case of the renormalization of both the electron and phonon 
propagation, i.e. by inclusion of the electron correlations n ktq , n k ^ q >, and N ktq<q /, we obtain 
the "self-energy" 

- illk,q(u) = — ~~~ ~ 227 (35) 

for the "propagator" D k>q (uj) in Eq. (0), where we approximated the electron self-energies 
by a constant damping factor. But this termination of the "continued fraction expansion" 
for the propagator D k ^ q (uj) is not the best one. Better results are obtained if the damping 
of the propagation of the scattered electron in n k) g,n kjqtq >, and N k ^ q i is included: 

- iU k>q {u) = -s — i2>y . (36) 

Now J2k'{ n k'-g ~ n k /)(t ) = for < q < q c , so that the propagator D kq is given as 
described in Ref. H by 



D k , q (t- t') = e(t-t') I 



A + e (-iA+-2 7 )(t-t>) X - e (-i\--2j)(t-t')\ ^ 




where 

1 f \ Vpq - Uq )2 + Ag 2^ . (38) 

The splitting of the phonon frequency into two modes as well as the above mentioned 
cancellation of certain corrections for the electron-hole bubble are special effects of the linear 
electronic energy dispersion. The terms proportional to k k q q r, K k qA ^ lk,q,q', and L k ^ q i in the 
time evolution of the phonon assisted density matrices h k>q and H k q (Eqs. fllfiD-flPTD), which 
led to a damping for h k>q and H k q) do not contribute to the time evolution of N q (Eq. QT5|)) 
due to the summation of h k , q over k. (There are contributions from electron states near 
the momentum cut-offs which do not cancel exactly. But in the spirit discussed above they 
may be neglected.) This cancellation occurs only due to the linear energy dispersion which 
makes the non-interacting energies e k — e k ~ q — u q = v Fq — ui q of h k , q and H k>q independent 
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of k. Therefore, the damping of h k , q and H kjq is only "existent" for their appearance in the 
differential equation (14) for n k , but not in that for N q . Although for the parameters chosen 
this remark is of little importance, it should have been taken into account for the equations 
in Born approximation with damping. 

The same kind of exact cancellation is responsible for the damping of only "half 
the electrons" in n k ^ q , n k ^ q >, and N k}qA i. These correlations enter the differen- 
tial equations for correlations of lower order only in the form n kA (i) := n kjq (t) + 

E q '>o( n k,q,q'(t)+n k+ql ^ ql (t)) = Efc' {{^l^l-q^k'-qYt + (Vi-g ^k ^k-q^k' ) % } and N k>q (t) '.= 
-E q > >0 (Nk, q ,q'(t) - N k+q ,^ q ,(t)) = iJ2k> {[^Wk-q^k'-qYt ~ {^l -q^k-q^k' )t] • ^ the 

non-interacting energy is the same for all these correlations, i.e. vanishes exactly due to 
the linear energy dispersion where momentum and energy conservation are identical for the 
electrons, only the quantities h kq and N kA must be considered. This fact was also used in 
the numerical calculations since it considerably reduces the number of differential equations 
to be solved. It also allows the extension of the transport equations by "half an order" , i.e. in 
the next order we take into account only certain sums of the correlations over momenta such 
that the number of indices does not increase. To be specific, we consider the correlations 
which arise from K k>q)q r, k k>q>q >, L k>q>q >, and lk, q , q ' by the substitution of one or both phonon 
operators B q by the density operator p q = J2k' Wk'-q^k'- These correlations as well as the 
extended transport equations are given in the Appendix. 

Figs. p| - |T4| show the numerical results for the Tomonaga-Luttinger model. As discussed 
in Ref. || the renormalization of the phonon propagation (Eq. fl37D ) is important for the 
Born approximation to describe the relaxation process adequately. The time evolution of the 
electron occupation number at the first phonon replica shown in Fig. [Tl| demonstrates that 
the kinetic equations including two-particle correlations (long dashed line) incorporate this 
renormalization in the dynamics as discussed above. For short times these equations give 
the same improvement for the time evolution of the second phonon replica in comparison 
with the Born approximation as in the one-electron problem (Fig. [12]). For larger times the 
agreement with the exact solution is worse because the important renormalization of the 
phonon propagation is not included in the differential equations for K k ^ q i, k k>q>q /, L kq)q >, and 
h,q,q'- This is achieved by the extension of the kinetic equations by the correlations in which 
the phonon operators B q are substituted by the electron density operators p q (Eqs. (|A1|)- 
( |A8|) ). The improvement by this dynamical renormalization can be seen in Fig. [12] where 
the results of the extended kinetic equations (named "highest order" , dot-dashed line) are 
in good agreement with the exact solution. Figs. [13] and [14| show the electron distribution 
function in the range of the first and second phonon replica for two different times. Observe 
the excellent agreement between the results of the kinetic equations of "highest order" (long 
dashed line) and the exact solution (solid line) even for large times (Fig. 111). 



C. Initial correlations 

Finally, we discuss the role of initial correlations for the dynamics. The initial correlations 
enter as initial values of the kinetic equations. In order to study the capability of the 
kinetic equations (|I~4"|)-(p4"|) to describe the effect of initial correlations we consider different 
initial states p inducing identical initial distribution functions n k (t ) and N q (t ). Since the 



9 



dynamics is exactly soluble for p = \ki, . . . , fcjv)(A;i, . . . , fcjv| the initial states were taken as 
convex combinations of states of this form. The exact solution is described elsewhere [13 . 
In the example shown in Fig. ^ the two initial states considered are 

(a) po = - {\k ,k - q B ,k - 2q B ,k - 3q B )(. . . | 

1 , 3 , 5 , 7 w ,1 
+|«o - gffB, k - ~q B , k - -q B , k - t^bX- • • I) 

( b ) Po = g (l^o, &q - 2<lB,k -q B ,k - ^fcX- ■ ■ I 

5 7 1 

+ - 2g B , fc - -q B , k - 3q B , k - -^Qb){- • • | } 

with a very small "resonant" phonon momentum q B = 2 ■ ^ in order to make the effect 
more pronounced. As the electron and phonon distribution functions n^ito) and N q (t ) 
as well as the phonon assisted density matrices h k ^ q (t ) and H k q (t ) are the same for both 
cases the Born approximation yields identical results whereas the kinetic equations including 
two-particle correlations provide an excellent approximation for the exact time evolution 
(Fig. 



IV. SUMMARY 

We studied kinetic equations for electron-phonon systems beyond the Born approxima- 
tion by the inclusion of many-particle correlations. The dynamics is described by a system 
of differential equations for the distribution functions and single-time correlations. The re- 
sults from the approximate kinetic equations were compared with the exact solutions for 
one-dimensional electron-phonon models. In contrast to real-time Green's function tech- 
niques, which have definite advantages in the context of quantum field theories, equations 
for retarded/advanced Green's functions are not needed. 

For systems with intrinsic damping, like the one-dimensional electron-phonon models 
studied, already the inclusion of two-particle correlations leads to a significant improvement 
in comparison with the Born approximation. The renormalization of the electron and phonon 
propagation, which is important to reduce the occurrence of unphysical results in the Born 
approximation, is automatically incorporated in the dynamics. There is no need for the 
Kadanoff-Baym ansatz. In addition to the renormalization of the propagators, the improved 
kinetic equations also correctly describe coherent two-phonon processes. We studied in detail 
the one-electron polaron model and the Tomonaga-Luttinger model. 

Finally, the ability of the improved kinetic equations to describe the effect of initial 
correlations was shown by comparison with exact solutions for the many-electron case. 
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APPENDIX: 

In this Appendix we present the extension of the kinetic equations mentioned in Sec. |111B 



In addition to the two-particle correlations, we consider the following correlations 

Pk,q,q'(t) = i Y {{BWk'-q^k-q'^k-q^k'Yt ~ (Bg> 1pU q 1p k -q'^k'-q)t) ( A1 ) 

w 

Pk,q,q>(t) = Y ((Bl4l- q 4-q>^-Mt + (BrfUt-rfk-rfu -q)tj (A2) 
k' 

rk,q,q'(t) = i J2 {{ B l4Wk-q-q4k^k'-q) C t ~ ( J B g /^_ ? V , |^fe-?- ? ''0fc')t) ( A3 ) 
k> 

Rk,q,q>(t) = J2 {{Bl^Wk-q-q'^k'-qYt + (B^-^IVfc-wV'k')*) , (A4) 
k' 

which arise from K k ^ qq ^ L kq ^/, and lk,q,q' by the substitution of one phonon operator 

B q by the density operator p q = J2k' Wk'-q^k', an d 



Sk, q>q >(t) = ((^k''i ) k'-qi ) k-q>^k-q^k'^k''-q')t + (^I'-q'^k'^l-q^k-q'^k'-q^k'^t) ( A5 ) 
k',k" 

Sk,q,q'(t) = % Y ((^l^l_ q ^l_ q 4k-q^k'^k''-q') C t ~ (i>k» -q'^k' ^k- q *P k-q'^Pk' -g^P k" )\ ) (A6) 
k',k" 

tk,q, q >(t) = ((^k''i ) k'i ) l-q-q>^k'-q^k''-q')t + (^I'-q'^k'-q^k-q-q'^k'^k^f) ( A7 ) 
k',k" 

W^'E ((i ) k''^k'i ) l-q-q''4 ) k^k'-q^k''-q')t ~ (i ) k''-q'^k'-q^li ) k-q~q'^k'^k'')f) , ( A8 ) 
k',k" 

in which the second phonon operator is also substituted by the corresponding density oper- 
ator. 

The extension of our transport equations by these correlations reads 

— fl k:q = . . . + Y 9q> {Pk,q,q> ~ Pk+q' ,q,q' + ^k,q,q' ~ ^k+q' ,q,q'} ( A 9) 
at q'>0 

~7,^k,q — ■ ■ ■ + Yj 9q' {Pk,q,q' ~ Pk+q' ,q,q' ~ Rk,q,q' + Rk+q' ,q,q'} ( A 10) 
dt q'>0 

{it + kk,Q,q ' = ' ' ' ~ 9 l' Pk >l''1 ~ 9qPk,q,q' ( A11 ) 

{it + 2 ^ Kk ' q,q ' = • • • + 9q' P k,q',q ~ 9qPk,q,q> ( A 12) 

{it + 2 ^ lk,q,q ' = " " ' + fl V rfc '9',<? + 9qn,q,q' ( A 13) 

{it + 2 ^ Lk,q,q ' = ' ' ' ~ 9q ' Rk > q '> q ~ 9qRk,q,q> ( A 14) 

where the dots stand for the terms in these equations written down so far, 
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(c^ + 2l ) Pk,q,q ' = ~^ Uq ' ~ VFq " >Pk ' q ' q ' 

+g q > {- (1 + Ng> - n k - q >) n k , q + (N q , + n fc _ g ) n k -^, q } 

+9q7i | { h k-q,qi ~ h k>q >) ^ h k' ,q ~ (H k -q,q> ~ H k>q >) ^ E k ^ q \ 
Z I k> fc' ) 

+9q^i n k'-q — nk')k k , q , q ' — g q 'S ktq>q i 



(A15) 



Jt + 2 V Pk >™ 



Jt + 2 V rk ' q ' q ' 



' d_ 

dt 



+ 27 R k , q , q 



(uy - v F q')p Kq , q > (A16) 
+9q' {- (! + N q' ~ n k - q >) N kt9 + (N q , + n k _ q ) N k _ q ,^ 

+fi'<?9 I (hk-q,q' - h k ,q>) Y Hk',q + (H k ^ q ^ - H k<q/ ) ^ ^k',q \ 
A I fc' fc' J 

+fi f <? ~~ n k')K kA>q ' + g q 'Sk,q,q' 

fc' 

-(uy - V F q')Rk,q,q' (A17) 

+fiV {- (! + N?' - ™fc- ? -<?') ™fc,<? + (^V + Wfc) ™fc-?',J 

+5 , '?9 I i h k-q,q' ~ ^fc,g') X! ^'-9 + ( H k-q,q> - #fc,g') ^ H k ^ q 1 
A I fc' fc' J 

~~ #9 ~~ n k')h,q,q' ~ 9q'tk,q,q' 

k> 

(uy - v F q')r Km , (A18) 
+5y {(1 + Ay - n k - q -f) N k:Q - (Ay + n k ) N k -q> iq } 

+9q^ \ -(h k -q,q> - h ktq >) ^ H k , )Q + {H k ^ q , - iffc, g ') ^ h k>,q \ 
Z I fc' fc' ) 

+9q XX nfe '-9 ~~ n k')L ky q yq ' + g q 'T k ^ q i 
k' 

(A19) 



and last but not least 



+ 



+ 27 ) Sfc ' 9 ' 9 ' = 9q ^^^ nk '- q ~ nfc ')Pfe,<?',<? (A20) 

{N k _ q ^ - N kjq r) H k',q + {n k - q ,q> - n k tf) h k',q \ + (<7 <~* <?') 

k' fc' J J 

+ 27^J S k ,q,q> = g q ^J2( U k'-q ~ ™fc')^fc,<?',<? (A21) 

- (q <-> g') 



1 

+ 2 



(A"fc_g,g' — A" fcj g/) /i fc / >? — [n k - q , q i — n ktq >) Y Hk',q 
fc' fc' 
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-{Nk-qrf — N ks i) H k / q + (n k _ q ^ q > — h ks i) ^2 hk',q 

k' k' 

j t + 27^) ?W = g q j- - n k >)R kt q',q 



1 

+ 2 



(N k -g y qi — Nk,g>) /ifc',g + (nk-q,qi ~ Tlk,q>) -fffc',q 



(A22) 
} + (?<- 

(A23) 



The same discussion as for the damping of the phonon assisted density matrices 
(Eqs. (p7D-fl34D) shows that the terms p..., P ., r .., and i? in the time evolution of and 



A^fc i9 lead to a damping of n k , q and iV^g by the same factor 27, as mentioned in Sec. [Ill A . 
In general, in all correlations containing a density operator p 9 , i.e. Y!,w(tyki-q$k' ■ • -)t> these 
electron propagations are not damped. 
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FIGURES 



k + q 



k + q 



FIG. 1. The interaction vertices for the polaron model. 





FIG. 2. The diagrams for the electron and the phonon distribution function. 




FIG. 3. The graphs for the time evolution of h kq and H k 



FIG. 4. The graph for the correlations K k ^ q /, k k ^ q >, L k ^ q >, and l k ^ q > 
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FIG. 5. The graph for the correlations nk, q , ^fc,q,q'j and N^^^t 



FIG. 6. Second order self-energy for the electron propagator. 
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FIG. 7. Time evolution of the electron occupation at the first phonon replica (k = ko — qs) 
for the one-electron polaron system with v = 1/16, T = 0.005, and q c /qB = 5. The exact so- 
lution is compared with the results from the kinetic equations in Born approximation with and 
without damping of the electron propagation and the improved equations including two-particle 
correlations. 
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The occupation at the second phonon replica (k=k -2q B ) 



n k (t) 




FIG. 8. Time evolution of the electron occupation at the second phonon replica (k = /en — 2qs) 
for the system of Fig. |?]. In addition to Fig. fj] the improved equations with damping are considered. 



The electron distribution function at co n t=20. 
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FIG. 9. The electron distribution function at uj^t = 20 for the system of Fig. f?|. Shown are the 
first and second phonon replica. The electron was initially at /en = 
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The electron distribution function at coJ=40. 



n k (t) 
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FIG. 10. The same as in Fig. 0, but for the time oj^t = 40. 
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0.2 



n k (t) 



0.1 



exact 

Born w. damping 

Born w. renorm. phonon 

w. corr. (damped) 




100 



FIG. 11. Time evolution of the electron occupation at the first phonon replica (k = ko — qs) for 
the Tomonaga-Luttinger model with v = 1/20, V = 0.003, and q c /qB = 3. The excited electron was 
initially at ko = 6qs {^f = 0). The exact solution is compared with the results from the kinetic 
equations in Born approximation with and without renormalization of the phonon propagation and 
the improved equations including two-particle correlations. 
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The occupation at the second phonon replica (k=k -2q B ) 
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FIG. 12. Time evolution of the electron occupation at the second phonon replica (k = — 2qs) 
for the system of Fig. |ll]. In addition to Fig. 11 the kinetic equations are extended by "half an 
order" . 



The electron distribution function at coJ=35. 
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FIG. 13. The electron distribution function at co^t = 35 for the system of Fig. 
the first and second phonon replica. The electron was initially at k$ = 6qB- 
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Shown are 
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The electron distribution function at coJ=100. 
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FIG. 14. The same as in Fig. [l3], but for the time uj^t = 100. The Born approximation without 
renormalization of the phonon propagation is omitted as it does not yield acceptable results at 



larger times (see Fig. 11 



The electron distribution function at co t=5. 




FIG. 15. The electron distribution function at u)$t = 5 for a four-electron system with v = 1/2, 
r = 0.01, and q c /qB = 10. The system was in two different initial states (a) and (b) described 
in the text with the same initial distribution functions. At time t = the electron distribution 
function was 0.5 at k = 33-36. 5qB and vanished elsewhere. The exact solution is compared with 
the results from the kinetic equations in Born approximation and from the improved equations 
including two-particle correlations. 
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